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Solid-state amorphization of crystalline copper nanolayers embedded in a Cu 64 Zr 36 metallic glass 
is studied by molecular dynamics simulations for different orientations of the crystalline layer. We 
show that solid-state amorphization is driven by a reduction of interface energy, which compensates 
the bulk excess energy of the amorphous nanolayer with respect to the crystalline phase up to a 
critical layer thickness. A simple thermodynamic model is derived, which describes the simulation 
results in terms of orientation dependent interface energies. Detailed analysis reveals the structure 
of the amorphous nanolayer and allows a comparison to a quenched copper melt, providing further 
insights into the origin of excess and interface energy. 
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I. INTRODUCTION 

Metals show a strong tendency to form crystalline 
phases and only appear in an amorphous state un¬ 
der certain conditions. Using far-from-equilibrium pro¬ 
cesses, metals can be kinetically trapped in a metastable 
state. This includes bulk metallic glasses (BMG), which 
are highly alloyed metallic systems quenched from the 
melt. The glass formation is supported by rapid quench¬ 
ing and the size difference of the component atoms. 1-3 
Furthermore, a crystalline metal sample can be forced 
into a disordered state by ion irradiation. The disor¬ 
der is introduced by high-energy impacts of ions, which 
disturb the ordered lattice due to local melt-quench 
processes. 4 Thin films produced with high deposition 
rates can also be amorphous. 5 In this case, the amor¬ 
phous state is metastable and only induced due to the 
high growth rate in the deposition process, while in equi¬ 
librium thin metal films on a variety of substrates usu¬ 
ally are crystalline. Examples include iron on amor¬ 
phous carbon substrates; 6 Cu, Ag, Al, Au, and Ni on 
sapphire substrates; 7-9 and Ni on tungsten substrates 10 
among many others. Amorphization is not limited to 
far-from-equilibrium processes, but can also happen for 
purely energetic reasons. High-angle, high-energy grain 
boundaries in some polycrystalline metal systems exhibit 
an amorphous structure due to the misorientation of the 
neighboring crystal lattices. 11 Molecular dynamics (MD) 
computer simulations on single-component 12 and binary 
alloy 13 Lennard-Jones systems identify a nanocrystalline 
instability. Nanocrystalline materials with grain sizes 
smaller than a critical value become unstable and collapse 
completely, leaving behind an amorphous metal, i.e., the 
grain boundary phase. Similarly, metallic nanoparticles 
below a critical size can also occur in an “amorphous” 
phase driven by the reduction of surface energy. 14 

In heterogeneous interfaces between crystalline metals, 
amorphous interphases were also found. 15 These inter¬ 


face interphases are thermodynamically stable and result 
from the misorientation and lattice mismatch between 
the adjacent crystallites. 16,17 This effect is called solid- 
state amorphization (SSA) and has recently also been 
discussed in the framework of complexion formation. 18 
Similar to the formation of interface interphases, a thin 
metallic film embedded in a different crystal phase can 
transform into an amorphous state if the thickness is be¬ 
low a critical value. 19 21 While energetically driven amor¬ 
phization of a thin crystalline layer due to size mismatch 
and misorientation to the abutting crystalline phases is 
a well known phenomenon, the amorphization of a thin 
elemental metal layer embedded in an amorphous matrix 
appears to be less likely, since the driving force should 
be significantly smaller. 

Ghafari et al ., however, showed recently that 
iron nanolayers embedded in an amorphous glass 
(Co75Fei2Bi3) can become amorphous, if the thickness is 
five monolayers (ML) or less, while at six or more mono- 
layers a crystalline phase is observed. 22 Whether or not 
this is a kinetic effect due the deposition conditions or an 
energetically driven phenomenon, which might depend on 
the lattice orientation, is the problem of interest in this 
study. 

In order to address this, we conducted MD simula¬ 
tions of Cu nanolayers of different orientation embedded 
in a Cu 64 Zr36 matrix and investigated the driving force 
behind the amorphization. Since there is no adequate 
iron alloy potential which also correctly models an al¬ 
loy that forms a metallic glass (MG), we instead used a 
system based on Cu and Zr. To investigate the thermo¬ 
dynamic stability of the amorphous layer, we start from 
crystalline nanolayers of varying thickness and check if 
there is a phase transition to the amorphous state at a 
critical thickness. We then develop a simple thermody¬ 
namic model based on the assumption that any amor¬ 
phization in this setup must be energetic in nature and 
check it against our simulation results. 
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FIG. 1. Time evolution of the composite Cu 64 Zr 36 -Cu systems with different nanolayer thickness. Exemplarily, we show a 
system in which a copper nanolayer with fee (100) surface was inserted. The top row (a) shows the evolution of the system with 
a copper nanolayer thickness of three monolayers. The nanolayer in this system amorphizes almost immediately. The middle 
row (b) depicts a system with four monolayers of copper, which stay partly crystalline. The bottom row (c) contains snapshots 
of a system where the nanolayer has a thickness of five monolayers. Here, the nanolayer stays crystalline. This simulation used 
the Mendelev potential. Copper atoms are shown in red, zirconium atoms in blue. Copper atoms that belong to the inserted 
layer are green, except for those in fee configuration, which are shown in yellow. 


II. SIMULATION METHODS AND SAMPLE 
PREPARATION 

A. Simulation method 

We conducted MD simulations using lammps. 23 The 
potential energy was modelled by a Finnis-Sinclair type 
potential by Mendelev et a/., 24 which was mostly fitted 
to the glassy state, and another potential by Ward et 
al. This second potential consists of elemental poten¬ 
tials by Zhou et al 25 that were combined into an alloy 
potential by fitting to several intermetallic phases using 
the Ward method. 26 To obtain independent confirmation 
of our simulation results, we ran simulations with both 
potentials. The time-step length was set to 2fs. In all 
simulations we employed periodic boundary conditions 
to obtain a surface-free system. 


B. Sample preparation 

Cu 64 Zr 36 metallic glasses were prepared by quenching 
from the melt at 2000 K to 300 K with a cooling rate 
of 0.01 K/ps. This procedure yields a MG with a local 
topology matching the experiment. 27,28 Two glasses with 
63108 atoms were prepared, one using the Mendelev po¬ 
tential and one using the Ward potential. The final size 
of the simulation box is approximately 10 x 10 x 10 nm 3 . 
All following steps were carried out twice, once with the 
Mendelev and once with the Ward potential. 

Copper nanolayers with the appropriate lattice con¬ 
stants at 300 K were created with (100), (110), and (111) 
surfaces. For each of these, the thickness varies between 
2 ML and 15 ML. To avoid stresses in the nanolayers af¬ 
ter insertion into the glass matrix, we scaled the x and y 
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(c) Layers with initial (111) interface. 


FIG. 2. Radial distribution functions of the copper nanolayers in several composite systems compared with the reference 
systems. Results of the simulation carried out using the Mendelev potential. 


dimensions of the glass to fit the nanolayers exactly and 
relaxed it again with a barostat applied in z direction at 
ambient pressure for 1 ns. 

The glass was cut at an arbitrary xy plane and the 
copper layers were inserted, so that the minimum ini¬ 
tial distance between any nanolayer atom and any matrix 
atom was at least 1.5 A. The possible consequences of low 
atomic distance and the resulting high potential energy 
at the interface are discussed later in section IV A. To de¬ 
velop a stable interface, the systems were equilibrated for 
1 ns at 300 K, again with a barostat applied only in the 
z direction. We kept the lateral dimensions constant be¬ 
cause any change in them would be dominated by the re¬ 
laxation of the MG. This would induce unwanted stresses 
in the nanolayer. At the end of this procedure, the sys¬ 
tems were completely equilibrated (cf. 29, section I). 

Further, we created reference systems of crystalline 
and amorphous copper phases. For the crystalline phase, 
we simply equilibrated fee copper at 300 K to obtain the 
correct lattice constant. Bulk amorphous copper was ob¬ 
tained by quenching from the melt at 2000 K with very 
high cooling rates. We had to employ a cooling rate of 
lK/ps for the Mendelev potential and 25K/ps for the 
Ward potential. These cooling rates are approximately 
the minimum cooling rates needed to avoid crystalliza¬ 
tion. The difference is a result of the different glass¬ 
forming ability of elemental copper in the two potentials. 

C. Analysis 

We applied a common neighbor analysis (CNA ) 30,31 as 
implemented in OVITO 32 to identify the structure of the 
nanolayers in the composite. The CNA calculates the co¬ 
ordination of all atoms by examining their neighborhood. 
To confirm these results, we calculated a radial distribu¬ 
tion function (RDF) by averaging RDFs calculated for 50 
snapshots of the equilibrated systems. The RDFs were 
calculated only for the atoms in the nanolayer. These 
RDFs were then compared to reference RDFs of the cop¬ 


per fee and copper glass systems. The short-range or¬ 
der of the amorphous nanolayers was analyzed using the 
Voronoi tessellation method implemented in OVITO, 32 
which divides the simulation cell into one polyhedron 
around each atom . 33-36 The polyhedra are characterized 
by the Voronoi index ( 723 , 714 , 715 , ne), where rii denotes 
the number of i-edged faces of the polyhedron. 

III. RESULTS 
A. Mendelev potential 

Figure 1 shows the time evolution of composite systems 
with different nanolayer thickness and the results of the 
CNA. Atoms depicted in yellow are nanolayer atoms that 
are fee coordinated. On insertion, the nanolayer had a 
(100) surface orientation. We can see that the nanolayer 
with 3 ML of copper amorphizes after a short simulation 
time, while the copper with 5 ML thickness stays crys¬ 
talline. At a thickness of 4 ML a mixed state occurs. 
The systems with different initial orientations show sim¬ 
ilar behavior. 

For obtaining an independent confirmation of the ap¬ 
pearance of an amorphous and crystalline copper phase 
at different nanolayer thicknesses, we calculated the 
RDFs of the nanolayers without including the glass ma¬ 
trix. The results are shown in Figure 2 and compared 
with the RDFs of the bulk crystalline and amorphous 
copper reference phases. Only three layers are shown for 
every initial surface orientation: the thickest amorphous 
layer, the thinnest crystalline layer, and the layer with 
a mixed state. The results match the CNA. Amorphous 
layers show a similar RDF to the bulk amorphous copper 
phase, including the characteristic double peak between 
4 A and 5 A. The RDFs of completely crystalline layers 
match the bulk fee copper, although the signal at large 
r gets smaller due to the finite size of the nanolayers. 
A mixed state is indicated by appearance of the second 
crystalline peak with reduced intensity. Depending on 
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(c) Layers with initial (111) interface. 


FIG. 3. Radial distribution functions of the copper nanolayers in several sandwich systems compared with the reference systems. 
Results of the simulation carried out using the Ward potential. 


the fraction of the amorphous phase, the glass double 
peak starts separating into the two clearly distinct crys¬ 
talline peaks (number three and four). In the system 
with initial (110) surface and a layer thickness of 5 ML, 
a small trace of crystalline phase is still detectable as in¬ 
dicated by a slight increase of the RDF at about 3.5 A, 
the position of the second crystalline peak. 

B. Ward potential 

The simulations using the Ward potential show the 
same behavior as the simulations using the Mendelev po¬ 
tential. The corresponding RDFs of the copper nanolay¬ 
ers are shown in Figure 3. The RDFs exhibit the same 
characteristics, i.e., the gradual disappearance of the sec¬ 
ond crystalline peak as well as the appearance of the glass 
double peak with decreasing thickness. The amorphiza- 
tion also occurs in slightly thicker layers. For the system 
with an initial (100) surface, 4 ML are amorphous, and 
even the 5 ML system only shows a small amount of crys¬ 
tallization. The mixed state for the initial surfaces (110) 
and (111) appears at 8 ML and 4 ML, respectively. The 
CNA results agree with the RDFs and are comparable to 
Figure 1. They are therefore omitted here. 

The simulations with both potentials agree that thin, 
initially crystalline nanolayers of copper become amor¬ 
phous. Thus, we can already exclude that amorphization 
is a result of kinetics in the deposition process, as the 
layers were inserted in a crystalline state. Still, the cre¬ 
ation of the interface may be connected with local heat¬ 
ing, leading to a fast melt-quench process in the smaller 
layers. It is therefore necessary to investigate the ther¬ 
modynamics of the system. 

IV. MODEL 

An indication for an energetically driven amorphiza¬ 
tion is already given by the fact that even initially crys¬ 


talline nanolayers undergo a phase transformation to the 
amorphous state. To explain the amorphization of the 
nanolayers and to test the hypothesis of energetically 
driven SSA, we propose a simple thermodynamic model. 
We formulate the internal energy U of the composite sys¬ 
tems. In the given ensemble, the free energy would be 
the appropriate thermodynamic potential, but given that 
the entropy term should favor the amorphization, we do 
not artificially increase the driving force for the amor¬ 
phization, but rather underestimate it. The internal en¬ 
ergy of a composite system with an embedded crystalline 
nanolayer is then 

U c = Umg + Nuf cc + 2A / y c -g. ( 1 ) 

Umg is the total internal energy of the bulk Cu 64 Zr 36 
glass phase, N is the number of atoms in the nanolayer, 
iq cc is the internal energy per atom of the copper fee 
crystal. Additionally, there are two interfaces, which con¬ 
tribute an energy of A. 7 c _ g each, where A is the interface 
area. 

If the system instead contains an embedded amorphous 
nanolayer, the internal energy is expressed as 

U a = Umg T T 2Ayg_g, (2) 

where u am is the per-atom internal energy of the glassy 
nanolayer and Ay g _ g is the glass-glass interface energy. 

Generally, it is to be expected that the internal energy 
of a copper crystal is lower than the internal energy of 
amorphous copper. As we nevertheless see amorphiza¬ 
tion of pure-metal nanolayers, the reason must lie in the 
interface energy. We propose that the crystal-glass in¬ 
terface energy 7 c _ g is higher than the glass-glass inter¬ 
face energy 7 g - g . In that case an amorphous nanolayer 
must be energetically favorable if its thickness doesn’t 
exceed a critical value, given that differences in interface 
energy can compensate the excess energy of the copper 
glass phase. A quantitative measure is provided by the 
internal energy difference, here expressed as an intrinsic 
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FIG. 4. Internal energies of the multilayer systems modelled with the Mendelev potential, (a) shows the internal energies 
extracted from the MD simulations as symbols and the linear regression results as lines. For visualization purposes a function 
h — h / • N/A with h! — 3eV was subtracted to exaggerate the difference in slopes between U c and U a . In (b)-(d), A U is plotted 
as a function of the number of monolayers for initial surface orientations (100), (110), and (111). Additionally, snapshots of the 
MD simulation are added. Here, red and blue atoms are copper and zirconium, respectively. Green atoms are copper atoms 
belonging to the nanolayer. Yellow atoms are fee coordinated. 


quantity independent of the surface area: 

A U = (U c - U a )/A 

N „ \ (N n \ 

^“^fcc + 27c-g ) — I ^-Wam + 27 g-g I ^2) 

= ^A^cu + 2A7. 

We can see that a negative value of AC/ signifies a sta¬ 
ble crystalline nanolayer, while a positive value of A U 


signifies a stable amorphous nanolayer: 

A U(N) < 0 crystalline nanolayer 

A U(N) > 0 amorphous nanolayer. 

Should the theory hold, we should be able to show that 
the critical number of atoms A^ cr it at AU(N cr [ t ) = 0 is 
the same as observed in the simulation by CNA. To get 
a more descriptive quantity, the number of atoms can 
easily be converted to the number of monolayers timl ot 
the thickness of the nanolayer d, as these quantities are 
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FIG. 5. Cut through a nanolayer of 3 ML thickness with an initial (111) surface at different time steps. Even after a simulation 
time of 100 ns, the mixed crystalline/amorphous state stays stable. The simulation was carried out using the Mendelev potential. 
Copper and zirconium atoms are shown in red and blue, respectively. Copper atoms belonging to the nanolayer are yellow if 
they are fee coordinated and green otherwise. 



FIG. 6. The nanolayer shown in Figure 5 removed from the metallic glass matrix and put in vacuum. The cut shows that the 
layer crystallizes almost immediately, proving that only the glass-glass interface stabilizes the amorphous phase. Yellow atoms 
are fee coordinated. 


proportional: 

N , 

— oc 71ml OC d. 

The missing parameters in our model are now iq cc , 
Wam, 7c-g, and 7 g_ g . Using equations 1 and 2, the internal 
energies of the different layer phases can be obtained from 
the relations 


nanolayer. This allows to calculate the interface energies 
either by directly using equations 1 and 2, or by subtract¬ 
ing Umg from the U-axis intercept of the U(N) curves. 
Both should yield the same result. We assume here that 
the interface energy is constant, see 29, section II for 
proof. 

A. Mendelev potential 


_dU c 

u fcc — i AT or u am 

d N 


dUa 
d N : 


(4) 


respectively. Alternatively, it would be conceivable to 
just use the internal energies of the bulk copper refer¬ 
ence systems. The problem would be that the amor¬ 
phous phase in the nanolayer is not necessarily the same 
as in the bulk. The bulk amorphous copper is quenched 
with very high cooling rates and has therefore more sim¬ 
ilarity to the melt. The amorphous copper phase in 
the nanolayer may exhibit different short-range order as 
it is allowed to relax to a low-energy state. Further¬ 
more, the ratio of interface to volume is very high, which 
means that the nanolayer structure is highly influenced 
by interface contributions. To calculate the interface en¬ 
ergy, the internal energy Umg is first taken from the 
pure Cu 64 Zr 3 6 glass sample before embedding the copper 


Figure 4a shows the internal energies U c and U a as a 
function of the number of atoms. All values are normal¬ 
ized to the interface area and Umg is already subtracted. 
The symbols show the internal energies extracted from 
the MD simulation, while the lines show the linear re¬ 
gression. The numerical data is listed in Table I. We 
note that all crystalline nanolayers have (approximately) 
the same interface energy in the Mendelev potential. The 
graph shows that the glass-glass interface energy is lower 
than the crystal-glass interface (in fact it is zero, which 
will be discussed in detail in section V), which fits the 
assumptions of our model. This lowered 7 g _ g favors a 
glassy nanolayer up until approximately 0.6atoms/A 2 , 
which corresponds to a critical thickness of about 

^crit ~ 7.5 A. (5) 
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FIG. 7. Internal energies of the multilayer systems modelled with the Ward potential, (a) shows the internal energies extracted 
from the MD simulations as symbols and the linear regression results as lines. For visualization purposes a function h — h! -N/A 
with h! — 3.35 eV was subtracted to exaggerate the difference in slopes between U c and U a . In (b)-(d), A U is plotted as a 
function of the number of monolayers for initial surface orientations (100), (110), and (111). Additionally, snapshots of the 
MD simulation are added. Here, red and blue atoms are copper and zirconium, respectively. Green atoms are copper atoms 
belonging to the nanolayer. Yellow atoms are fee coordinated. 


The thickness d cr [ t can only be given approximately, due 
to the different densities of the two phases and the rough 
interface. Figures 4b-d show AC/ as a function of the 
number of monolayers. The direct comparison reveals 
that the calculated energy differences and the solid-state 
amorphization observed by the CNA method agree very 
well. At the critical thickness A U = 0 we observe a mixed 
crystalline/amorphous nanolayer. The figures show that 
the transition is not as sharp as our model assumes, a 
partly crystalline layer also exists for A V slightly greater 


than zero. That is a result of omitting a description 
of the two-phase region in the model; entropy and the 
additional interfaces play a role here. Nonetheless, the 
critical thickness is correctly reproduced without these 
complications. A further comparison with the RDFs in 
Figure 2 leads to the same conclusions. While the good 
fit of the model with simulation data supports the con¬ 
clusion that the SSA is due to energetic reasons, we also 
investigated the influence of the simulation setup on the 
amorphization of the nanolayers in 29, section III. We 






































find that the interface creation leads to a heat spike, but 
that this only serves as activation energy for the SSA 
process. Crystalline layers with a thickness below d cr it 
can be produced, but are not energetically favorable. 

In order to verify that the amorphous phase is indeed 
stable over a long time scale and that this is a result 
of the glass-glass interface, we conducted two further 
simulations. In the first simulation, we simply took a 
composite system with a mixed crystalline/amorphous 
state in the copper nanolayer and let the simulation run 
for 100 ns. If the amorphous phase is only produced 
by, e.g., stress in the initial system after insertion of 
the nanolayer, the crystalline phase should start grow¬ 
ing again over the longer time frame. The simulation 
results are depicted in Figure 5 and show that the mixed 
state is stable, as predicted by our model. 

A direct proof that the solid-state amorphization is due 
to the presence of a glass-glass interface was obtained 
by removing the glass matrix. The results are shown in 
Figure 6. The free amorphous layer crystallizes almost 
immediately, as would be expected. 


B. Ward potential 

Figure 7a shows the internal energies as calculated with 
the Ward potential as a function of the number of atoms. 
All values are again normalized to the interface area and 
Umg is already subtracted. The symbols show the inter¬ 
nal energies extracted from the MD simulation, while the 
lines show the linear regression. In the Ward potential 
the (111) interface has a slightly lower interface energy 
than the (100) and (110) interfaces, which are approxi¬ 
mately the same (see also Table I). As with the Mendelev 
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FIG. 8. Energy-volume curves for amorphous and crys¬ 
talline copper and data points for amorphous and crystalline 
nanolayers. Atomic volumes of the nanolayers were obtained 
by a Gauss fit to the Voronoi volume distribution. The crys¬ 
talline nanolayers vary in density with their thickness, ap¬ 
proaching the equilibrium volume of an fee crystal at higher 
thickness. 


potential, the glass-glass interface energy is lower than 
the crystal-glass interface energy, favoring an amorphous 
nanolayer up to the critical thicknesses 


j(100) , 

u crit " 

« 10.8 A, 

(6) 

d (110) f 

“crit 

» 10.7 A, and 

(7) 

d (111) F 

“crit 

a 9.9 A. 

(8) 


The difference in critical thickness is a result of differ¬ 
ent 7 c _g for the three surface orientations. The transi¬ 
tion thickness is higher than in the Mendelev potential 
despite a smaller Ay, as the excess energy of the amor¬ 
phous phase is lower. By plotting AU as a function of the 
number of monolayers, a direct comparison to CNA and 
RDF results is possible. Figures 7b-d show AU(timl) 
compared with snapshots from the simulation. Again, a 
good match between the nanolayer phases shown in the 
snapshots and the predicted critical thickness is visible. 
For the same reasons as stated earlier, a mixed state oc¬ 
curs. 

All in all, the results using both potentials agree qual¬ 
itatively and support our thermodynamic model. There¬ 
fore, a purely kinetic reason for the amorphous nanolay¬ 
ers can be ruled out and an energetic picture of solid-state 
amorphization can be supported. 


V. STRUCTURE AND ENERGY 
A. Mendelev 

We compared the energy of the amorphous nanolayer 
copper phase with the reference bulk amorphous cop¬ 
per phase. The energy-volume curve in Figure 8 shows 
that the nanolayer phase is energetically higher than the 
bulk phase. Additionally, an expansion of the crystalline 
layer at low thicknesses is visible, possibly due to in¬ 
terface stress. Both the peculiarity of the zero glass- 
glass interface energy, as well as the higher excess energy 
of the amorphous phase can be linked to the structure 
of the phase. Figure 9a shows the Voronoi statistics of 
those copper atoms in the nanolayer, that are only sur¬ 
rounded by other copper atoms. This allows a compar¬ 
ison to bulk amorphous copper: in contrast to the bulk 
phase, the nanolayer phase contains no twelve-fold co¬ 
ordinated atoms. This higher energy structure is stabi¬ 
lized by the interface, as shown in Figure 9b: the Voronoi 
statistics of the whole nanolayer (including copper atoms 
that have zirconium neighbors) are very similar to the 
bulk Cu 64 Zr 36 . This reduces the interface energy to al¬ 
most zero. 
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FIG. 9. Voronoi analysis of the amorphous nanolayers in sys¬ 
tems simulated using the Mendelev potential. In (a) the bulk 
amorphous copper phase is compared with the amorphous 
copper phase in the nanolayers. Only those nanolayer atoms 
were included, which were surrounded solely by other copper 
atoms. In (b) all nanolayer atoms (including those with zirco¬ 
nium neighbors) were considered and compared to Cu 64 Zr 36 
bulk. 


FIG. 10. Voronoi analysis of the amorphous nanolayers in 
systems simulated using the Ward potential. In (a) the bulk 
amorphous copper phase is compared with the amorphous 
copper phase in the nanolayers. Only those nanolayer atoms 
were included, which were surrounded solely by other copper 
atoms. In (b) all nanolayer atoms (including those with zirco¬ 
nium neighbors) were considered and compared to Cu 64 Zr 36 
bulk. 


B. Ward 

In the Ward potential the amorphous nanolayer also 
features a structure with different energy than the bulk 
amorphous copper (—3.39 eV/atom for the reference sys¬ 
tem, —3.40eV/atom for the amorphous nanolayer). The 
explanation can again be found in the Voronoi statis¬ 
tics of nanolayer copper atoms surrounded completely by 
other copper atoms: Figure 10a shows that the twelve¬ 
fold coordinated atoms are missing again. The com¬ 
parison of the Voronoi statistics of the whole nanolayer, 
though, show a difference to the Cu 64 Zr 36 MG, especially 


concerning the (0,0,12,0), (0, 2, 8, 2), and (0, 3,6, 3) 
polyhedra (Figure 10b). This leads to a small but non¬ 
zero glass-glass interface energy. 

VI. CONCLUSIONS 

Using MD simulations, we observed the amorphization 
of elemental copper nanolayers embedded in a Cu 64 Zr 36 
metallic glass if the layer thickness stays below a critical 
value. This is in accordance with experimental results, 
which report thin amorphous iron nanolayers embedded 
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TABLE I. Internal energies and interface energies extracted from MD simulations and predicted critical thicknesses. 


Potential 

Initial 

interface 

^fcc 

'Warn 

Aixcu 

7c-g 

7g-g 

Ay 

N cr it / A 
(atoms/A 2 ) 

d C rit 

(A) 


(eV / atom) 



(meV/A 2 ) 
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(110) 
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0.846 

10.7 


(111) 

-3.46 

-3.40 

-0.06 

27.6 

5.5 

22.1 
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9.9 


in Cc> 75 Fei 2 Bi 3 . 22 We could show that the amorphization 
is not a kinetic effect due to deposition, as our simulations 
start from a crystalline state. Rather, the glass-glass in¬ 
terface energy is significantly lower than the crystal-glass 
interface energy, which stabilizes the amorphous copper 
phase. This solid-state amorphization is similar to the 
case at heterogeneous crystal interfaces, except that in 
our case the reduced glass-glass interface energy is suffi¬ 
cient to induce amorphization. At a critical layer thick¬ 
ness, which is on the order of a nanometer, a mixed crys¬ 
talline/amorphous state appears. This state is also stable 
over longer times, which further supports the picture of 
solid-state amorphization: if the amorphous state is only 
a result of stresses in the initial setup, the crystallites 
in the layer would grow again with time. They instead 
keep their size. Analysis of the amorphous structure in 
the nanolayer further confirms that the interface energy 
is a dominating factor in the structure of thin nanolay¬ 
ers: if it can be reduced, the amorphous layer can even 
be driven into a state with a higher bulk energy than 


a quenched melt. Technological applications for glass- 
glass composite systems have already been proposed in 
the realm of magnetic tunnel junctions 22,37 and could 
benefit from further research into different multilayer sys¬ 
tems. 
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Supplemental Material on the Simulation Setup 

Solid-state amorphization of Cu nanolayers embedded in a Cu 6 4Zr 36 glass 

This supplemental material discusses some additional technical aspects of the simulation setup in 
detail. In the first section, we show that the systems are equilibrated. We demonstrate in the 
second section that the interface energy is a constant independent of the layer thickness. In the 
final section, we discuss the kinetics during interface creation and show that they do not influence 
the thermodynamics of the system. 


I. EQUILIBRATION OF GLASS MATRIX AND 
COMPOSITES 

To extract the correct energies from the molecular dy¬ 
namics (MD) simulations, it is necessary to make sure 


that the systems are equilibrated, i.e., that the energy 
of the system is constant in time. The composites were 
equilibrated for 1 ns at 300 K, which is sufficient to equi¬ 
librate the composite systems. This is shown exemplarily 
for the Mendelev potential, for systems with initial (100) 
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FIG. 1. Internal energy over time for systems simulated with 
the Mendelev potential, (a-b) Internal energy of two metallic 
glass-copper composite systems, (c) Internal energy of the 
pure metallic glass. 


FIG. 2. Internal energy over time for systems simulated with 
the Ward potential, (a-b) Internal energy of two metallic 
glass-copper composite systems, (c) Internal energy of the 
pure metallic glass. 
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interface and 2 monolayers of copper (Figure la), as well 
as 15 monolayers of copper (Figure lb). For the Ward 
potential, the same is shown in Figures 2a and 2b. To 
ensure that the energy of the glass matrix stays constant 
over the simulation time, we ran simulations with only 
the glass matrix for 1 ns. These are shown in Figures lc 
and 2c and show that the glass is already equilibrated 
and does not change in internal energy. 


II. INDEPENDENCE OF THE INTERFACE 
ENERGY FROM LAYER THICKNESS 

Our model assumes an interface energy independent 
of the layer thickness. To prove that this assumption is 
correct, we additionally performed simulations on a sys¬ 
tem with 100 monolayers of copper using the Mendelev 
and the Ward potential. We inserted the results of this 
simulation into plots of the internal energy as a function 



N/A (atoms/A 2 ) 


FIG. 3. Internal energy over number of atoms for crystalline 
nanolayers with initial (100) surface in the Mendelev poten¬ 
tial. The blue circles show the results from Figure 4a in the 
main paper. The red circle is the system with 100 monolayers 
of copper. To exaggerate the slope for visualization purposes, 
we subtracted a function h = 3eV • N/A from U/A. 


TABLE I. Comparison of the interface energy in a composite 
system with initial (100) surface modelled using the Mendelev 
potential. 

Uf cc Umg/ Nm g 7c-g 

(eV/atom) (eV/atom) (meV/A 2 ) 

5-15 layers ^T20 ^440 36A 

100 layers —3.20 —4.40 35.6 


of the number of layer atoms, both normalized to the 
interface area. Figures 3 and 4 show that the compos¬ 
ites with very big copper layers have the same energy 
U c as predicted by extrapolation using our model (line). 
Furthermore, we calculated Uf cc for the copper nanolayer 
atoms far from the interface, as well as the per-atom en¬ 
ergy of the metallic glass matrix far from the interface 
and compared them to the values from Table I in the 
main paper. Tables I and II list the results and clearly 
show that the results in a system with “bulk” copper 
agree with the results for composites with thin copper 
layers. The small deviation of the interface energy, as 
well as the deviation of Umg for the Ward potential, he 
inside the range of numerical errors. The independence 
of the interface energy from the layer thickness is already 
a reasonable assumption, given the linear behavior of U c 
and U a . With these additional data, the assumption is 
proven. 



N/A (atoms/A 2 ) 


FIG. 4. Internal energy over number of atoms for crystalline 
nanolayers with initial (100) surface in the Ward potential. 
The blue circles show the results from Figure 7a in the main 
paper. The red circle is the system with 100 monolayers of 
copper. To exaggerate the slope for visualization purposes, 
we subtracted a function h = 3.35 eV • N/A from U/A. 


TABLE II. Comparison of the interface energy in a composite 
system with initial (100) surface modelled using the Ward 
potential. 

Ufcc Umg/ Nm g 7c-g 

(eV/atom) (eV/atom) (meV/A 2 ) 

6-15 layers ^3T6 ^448 29T 

100 layers —3.46 —4.49 29.3 
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FIG. 5. Temperature profiles for composite systems simulated using the Mendelev potential. The top row shows a system with 3 
monolayers, the bottom row a system with 7 monolayers of copper. The graphs are produced by averaging over 2 A thick slices, 
of which there are 500 in z direction. The top part of the graphs show the ratio of layer atoms to matrix atoms in the current 
slice to indicate the position of the two phases and their interface. The lower part of the graphs show the temperature in the 
current slice. Below the graphs, the corresponding snapshots of the composite systems are shown. The leftmost snapshots are 
always from the simulation step with the highest temperature spike. Note: After creating the composite systems, we initialize 
the momenta from a Maxwell-Boltzmann distribution. For the first few steps the distribution between kinetic and potential 
energy is not equilibrated, which leads to an instantaneous temperature higher than 300 K in the glass matrix. 


III. CREATION OF THE COMPOSITE 
SYSTEMS AND KINETICS AT THE INTERFACE 

After inserting the copper nanolayers, some atoms 
from the nanolayer may have a small distance to atoms 
from the metallic glass matrix. This could lead to a 
high potential energy and consequently initial temper¬ 
ature spikes at the interfaces, which may be responsi¬ 
ble for the amorphization of the nanolayers. In fact, for 
some time at the beginning of the simulation, tempera¬ 
ture spikes appear mainly on the matrix side of the in¬ 
terface (Figures 5 and 6). In case of the systems with 
nanolayers thick enough to stay crystalline, the temper¬ 


ature inside of the layer stays comparatively low, while 
there is an increase of temperature in the layers when 
they become amorphous. While this may be a sign of the 
amorphous phase being the result of a fast melt-quench 
process, it may also be that the initial heat spike only 
provides an activation energy for an energetically favor¬ 
able solid state amorphization. To test this, we used 
a modified simulation setup: The system was prepared 
in the same way as described in the main paper, but the 
nanolayers were held fixed in the crystalline state for 1 ns. 
This avoids penetration of heat into the layer. After that, 
the nanolayers were released and the whole system equili¬ 
brated normally for an additional 1 ns. Additionally, we 
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FIG. 6. Temperature profiles for composite systems simulated using the Ward potential. The top row shows a system with 3 
monolayers, the bottom row a system with 7 monolayers of copper. The graphs are produced by averaging over 2 A thick slices, 
of which there are 500 in z direction. The top part of the graphs show the ratio of layer atoms to matrix atoms in the current 
slice to indicate the position of the two phases and their interface. The lower part of the graphs show the temperature in the 
current slice. Below the graphs, the corresponding snapshots of the composite systems are shown. The leftmost snapshots are 
always from the simulation step with the highest temperature spike. 


ran molecular statics simulations on the initial system, 
which also avoids heating the system locally. The results 
are shown in Figures 7 and 8. The molecular statics sim¬ 
ulations never lead to an amorphization of the copper 
layer. The reason for that is simply that a phase trans¬ 
formation from fee copper to amorphous copper needs a 
certain activation energy. We can also see that the energy 
of the statically relaxed systems is always much higher 
than the MD-equilibrated systems. This shows that even 
the creation of a low-energy crystal-glass interface needs 
some activation energy. When holding the copper layers 
fixed, very small layers will still become amorphous, while 
the amorphization is suppressed in thicker layers. These 
crystalline nanolayers are metastable, as can be seen in 
the graphs: Their energy agrees with the prediction of 


our model (blue line), which also states that A U > 0 
in this case. Only above the critical thickness, at the 
crossover of U a (dashed line) and U c (blue line), are the 
crystalline nanolayers stable again. This is in complete 
agreement with our model and shows that initially fix¬ 
ing the copper layers only kinetically inhibits the phase 
transformation on the simulated time scales. 
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FIG. 7. Internal energy of the composite systems as a func¬ 
tion of the number of atoms in the nanolayer. This graph 
shows nanolayers with initial (100) surface in the Mendelev 
potential. Data was normalized to the interface area and a 
function h = 3eV • N/A was subtracted from U/A to exag¬ 
gerate the difference in slopes for visualization purposes. The 
blue circles show the results from Figure 4a in the main paper. 
The yellow diamonds show the results of an MD simulation 
where the layer atoms are first held in position, followed by a 
normal MD simulation. Except where marked, the nanolayers 
stayed crystalline in these simulations. Red squares show a 
simple molecular statics calculation plus a 2 x 3 / 2 NksT term, 
to compensate for the missing temperature contribution to 
the internal energy. 



FIG. 8. Internal energy of the composite systems as a function 
of the number of atoms in the nanolayer. This graph shows 
nanolayers with initial (100) surface in the Ward potential. 
Data was normalized to the interface area and a function 
h = 3.35 eV ■ N/A was subtracted from U/A to exaggerate 
the difference in slopes for visualization purposes. The blue 
circles show the results from Figure 7a in the main paper. 
The yellow diamonds show the results of an MD simulation 
where the layer atoms are first held in position, followed by a 
normal MD simulation. Except where marked, the nanolayers 
stayed crystalline in these simulations. Red squares show a 
simple molecular statics calculation plus a 2 x ^ANksT term, 
to compensate for the missing temperature contribution to 
the internal energy. 




